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Method of Moments for computing electromagnetic 
scattering from a rough cylinder 

Rahul Trivedi and Uday K Khankhoje 


Abstract 

In this tutorial paper, we formulate a two-dimensional integral-equation based method of moments approach for numerically 
computing the electromagnetic fields scattered from an azimuthally-rough dielectric cylinder or an axially-rough perfectly conduct¬ 
ing cylinder. The electric field integral equation is discretized in both the cases to yield a system of linear equations, and a detailed 
evaluation of the matrix elements involved in this system of equations is presented. Numerical computation of important quantities 
such as the scattering cross-section or the scattered electric fields is also discussed. The primary purpose of this numerical scheme 
is to validate perturbative solutions that arise in the context of rough cylinder scattering problems. 


I. Introduction 

Evaluating the electromagnetic scattering from cylinders with rough surfaces is a computationally tedious task. This is because 
the stochastic nature of the surface necessitates taking an ensemble average of the scattered fields over many realizations of a 
rough surface till suitable numerical convergence is achieved. In spite of their accuracy, full-wave techniques such as integral- 
equation Q or finite-element Q based methods, are impractical in remote-sensing applications where the scattering from a 
variety of cylinder configurations needs to be estimated Q, Q. Thus, it is natural to formulate perturbative solutions, which, 
for some tradeoff in accuracy, provide fast estimates of the scattering solution. 

A second order perturbative method for computing the electromagnetic fields scattered by a homogeneous cylinder with a 
rough surface has recently been proposed 0.@. In order to validate such a solution, it is essential to compare the perturbative 
solution with the solution predicted by a full-wave technique. In this paper, we describe the formulation of an integral-equation 
based method of moments which is used for the purpose of validating the above referenced perturbative solution. 

While the most general formulation of this problem is in three-dimensions, we restrict ourselves to two possible two- 
dimensional cylinder configurations - the extension to three dimensions is straightforward, but beyond the scope of the present 
paper. In the first configuration, we consider scattering from a dielectric cylinder of infinite length whose roughness function 
has only azimuthal dependence. By employing the two-dimensional Green’s function in a homogeneous medium, followed by 
an expansion of unknown currents in triangular basis functions and testing along pulse functions, we formulate our scattering 
solution in Section In the second configuration, treated in Section we consider scattering from a perfectly conducting 
cylinder of finite length with a roughness function which is only axially-dependent. Unlike the previous case, we require the 
three-dimensional Green’s function here, and formulate the solution by expanding the unknowns in pulse basis functions and 
test along delta functions. Detailed remarks on the discretization parameters and evaluation of the associated matrix elements 
is made in both cases. 


II. Scattering from an azimuthally rough dielectric cylinder 

Consider an azimuthally rough cylinder with it’s surface being described by r{(j),z) = rs{(l)) = a + h{cj)) (where {r,cj),z) 
are the cylindrical coordinates of the point in question), and as illustrated in the schematic in Fig. [T] The dielectric constant 
of the cylinder is assumed to be and it’s permeability ^d- It is of interest to calculate the scattered fields when a plane 
wave is incident on this cylinder. We consider only TM polarisation (E||z and H _L z), which corresponds to an incident field 
with a = 7r/2 and i ?2 = 0 in Fig. A harmonic time dependence of the form exp{jujt) is assumed in all the fields and 
suppressed throughout this tutorial. The incident field is given by: 

E™'^ _ _ exp(—jfco(a:cos(/)o + ysini^o))-^ = exp(—jfcor cos((^ — (j)o))z (1) 


which is a plane wave propagating at an angle (pg with the x axis and polarized in the 2 ; direction. The starting point of our 
analysis is the electric field integral equation (EFIE) for homogeneous dielectric scatterers ||^ Ch. 1.9]: 
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Fig. 1: Schematic of the rough cylindrical scatterer with a stochastic roughness h{(j>,z) and illuminated by a plane wave 
propagating with wavevector k; in the xz plane at an angle of a with the x axis. Ei is the component of the electric field in 
the xz plane, while E 2 is the component along the y axis. 


where = A^f'^z and = Fx^'^x + Fy^'^y for p G {0, d} are the magnetic and electric vector potential respectively 
defined by (the superscripts indicate the physical region in which the expression is valid; 0 corresponds to vacuum and d 
corresponds to dielectric): 

2-k 

^ (3a) 


A‘f \x,y) = ^ J + r'^{(!>') dcj)' 

0 

27T 

F^P\x,y) = dc^' 


(3b) 


where J{(p) and K{(j)) are the equivalent electric and magnetic currents along the cylinder interface, and R = [{x — 
rs{4>') cos 4>')‘^ + {y — is the distance between a jjoint {rs{4>'),(j)') on the cylinder surface to the point 

of observation {x,y). t{(j)) denotes a unit vector tangential to the cylinder surface at the point (rs{4‘),4>) nnd is given by: 




x{—rs{(l)) sin 0 + rg{(j)) cos(j)) + y{rs{4') cosc^ + r'g{(j)) sinc^) 




(4) 


Using Eqs. and Q, it can be shown that: 
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(5) 


To numerically implement a MoM routine to solve the EFIE (Eq. 0 )^ we assume a piecewise linear approximation of the 
cylinder surface (in 0) between N discrete points {rn+i/2,4>n+i/2), where 4>n+i/2 = + 5 ) and rn+ 1/2 = rsi'pn+ 1 / 2 ) 

V n S {0,1,2,... A — 1}. Thus, in this approximation, the radial distance of a point on the cylinder surface at an azimuthal 
angle (j) satisfying (j)n-i /2 < </>< </>„+i /2 is given by: 
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where, for notational convenience, we have defined r„ = [vn+xli + ^n- 1 / 2)/2 and = ^(r„+i /2 — rn-ii 2 )- In the analysis 
that follows, we denote the discretization angle I-k jN by 8. 

The unknown functions and are then approximated as a linear combination of triangular basis functions (shown 

in Fig. 1^; 


Af-l 


AT-l 


m « E 
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where Vn € {0,1,2 ... — 1}: 




1 — ^\(j> — n6\ for {n — 1)5 < (j) < {n + 1)S 


0 otherwise 

Further, to impose Eq. (|^, we employ the pulse testing functions Pn{4>) (see Fig. defined by (Vn £ {0,1,2,... At — 1}): 

1 for (n — 1)5 < 4> < (n + ^)5 


( 8 ) 


PnW = . 

I 0 otherwise 

Consider now the inner products of the individual terms in Eq. (j^ with Pm{4>)- For instance; 

27r 27r 27r 

I I JmPm{m^o\k^R)\l’rW) + r'i{<i2') d<5 d^' 

0 ^ 0 0 

jV—1 . N—1 

= E ^ / / ^n{(l)')Prn{4>)H^Q\kpR)^ d(j) dfj)' = E 

n=0 rn=0 


(9) 


where 


(m+l/2)S (n+l)5 


Q,(P) = 

^m.n 


4j 


^n{(l>)H'^\kpR)\jrl{4>') + r'^{(t)') d<j) dcj)' 


( 10 ) 


( 11 ) 


{m —1/2)5 (n—1)(5 


wherein we have used the fact that Pm{4>) is supported only on ((to — 1/2)5, (to + 1/2)5) and A„(/)) is supported on 
((n — 1)5, (n + 1)5). It is usually not possible to obtain a closed form expression for the integral in Eq. ( [IT) . However, for 
sufficiently small 5, the integral aZ/n can be approximated as: 


v(p) 
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Eig. 2; Triangular basis functions, A„((/)) and pulse testing function, Pn{4>) used in the MoM implementation. 
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with 


(m+l/2)i5 {n-l/2)5 
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(m+l/2)(5 (n+l/2)(5 


Ar,{c^')Hl/\kpR) d(^ d^' 
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^ / / 
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where r„_ 3/4 = ( 3 r„_i /2 + r„- 3 / 2 )/ 4 , r„+ 3/4 = ( 3 r„+i /2 + rn+ 3 / 2 )/ 4 . In writing Eq. (12i, we have split the integral 


(j)' into three different integrals and approximated yr'g{(j)') + r'g{(j)') centroidally in each of them. Depending upon the values 
of m and n, the integrals may or may not be singular. We deal with them in four cases: 

1) For m ^ n, n— 1, n +1, the integrands in Eq. ( [T^ are not singular, and for sufficiently small 5, they can be approximated 
centroidally: 
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Hq (^p-^m,n+3/4) 
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where = y'r^+r^- 2 r„r„ cos((m - n)S). 

2) For m = n, the integrands in Eq. ( [T3] l are singular and have to be integrated analytically. To do so, we use the small 
argument approximation to the Hankel functions: 


Hq^\x) « 1 — — In ( -^ ) for a; —)■ 0 


2j f jx 


(15) 


where 7 « 1.78107 is the Euler’s constant. Consider the integral in 2 - Using Eq. §’ 
Together with the small argument approximation, the integral becomes: 
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In \p — P'\Aq{P') dp dp' 


The integral in Eq. can be evaluated analytically to give: 
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A similar procedure can be adopted to evaluate i and 3 . We give the hnal result below: 
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I ln(S) - 0.0699 


2+1/2 4 

3) For m = n + 1, a^^-i 1 1*^ evaluated using the centroidal approximation, since the integrand is not singular: 
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(p) 

A procedure similar to that used in evaluating ^ 2 used to evaluate am, 771 - 1,2 and 
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4) The case of m = n — 1 is identical to the previous case; 
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Consider now the inner product of 



dF^P'l 1 

dx 
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with ((/>): 


" " PM 
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( 22 ) 


where 

(m+l/2)5 (n+l)(5 

/3m/ = ^ cos{(j) - (/>') + Pg(^^’’S(</>) sin((()' -(/>)- rs(«/')) dcj)' d(j) (23) 


(m—1/2)5 (n—1)5 


which we again express as a sum of three integrals: /3m/ = (/3// ^ + /3^/ 2 + /3m/ 3 ) where; 


(m+l/2)5 (n+l)5 


/ 3//.1 = / / A„(<^') 
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F{(j), (/>') dcj) 
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k^7n,7l,2 


R^p) 

r^7n,7l,3 


M^^^\kpR) 


d^' del 


(24a) 


(24b) 


(m—1/2)5 (n—1/2)5 
(m+l/2)5 (n—1/2)5 


( 2 ) ^ 

A» (<//') ^° del' del 


(24c) 


(m—1/2)5 (n —1)5 

where F{(j), (p') = rs{4>')rs{4‘) cos{(f) — ep') + r'g{(p')rs{(p) sin((/) — ep') — rl^ep'). We again consider four different cases: 
1) For m ^ n, n — 1, n + 1, the integrals in Eq. ([24)i are not singular and can be approximated centroidally: 
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r-^m,71,1 
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m,n+3/4 


ikMn+ 3 / 4 ) 
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Pm,71, 3 Q TD ‘ ^m, 71— 314 

o ^m,71—3/4 


(25) 
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where Fm,n = ^'m^n cos((n — m)J) + rmS„ sin((m — n)(3) — with s„ denoting the slope of the segment containing 
r„. For instance s„_ 3/4 = (r„_i /2 - r„_ 3 / 2 )/^- _ 


2) For m = n, the integral is singular. Additionally, note from Eq. (2 1 that what needs to be evaluated is not 


dx dy 
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rather [—g|-g 5 ^]ls+ [^g|-the integrand is not singular, then the limit is the same as the value 
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of the inner product at the surface S, however care must be taken in computing the limit for the case when the integrand 
is singular. Consider, for instance, the evaluation of 2' 

{m-\-l/2)S (m+l/2)(5 
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,(o) 

m,m,2 


= lim 


An(A) ^° d^' dcj, 


R 


(m—1/2)(5 (m—1/2)5 


with 


R = x/+ As(0) + e)^ - rs{(t)'){rs{4’) + e) cos((^ - (j)') « 


(28) 


(29) 


wherein we have used the fact that since \(f) — (j)'\ <5, \4' — (j)'\ is small and thus cos{(f> — 0') « 1 — ^ ^ and 

sin(^ — (j)') Ki {(f) — (f)'). Also note that e > 0, since 
use the small argument approximation and set (koR) « — ■ Moreover, for small |()) — (j)'\, the term 


arises in the evaluation of [—-g^]|5+. We may also 


can be approximated by [erm - ^{4' ~ 4’’Y]- Therefore: 
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-5/2 -S/2 

where we have performed the substitution (js — mS 4> and (j)' — mS 4' in th® step. Consider the integral: 
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Since, in Eq. (30i, \4'\ < <5/2, 


s /2 


lim 

£->■0 


r'L{4-4'Y - 2 e?'„ 


-5/2 


rl,{4-4’f 


•d4 = 5 — 2Tr 


(31) 


(32) 


Note the extra 27r term in the above integral. Had we taken the limit before evaluating the integral, we would have 
obtained <5 as the result of the above integration, which would yield erroneous results. Substituting this integral into Eq. 
(Boll, we obtain: 

,(o) _ 3j 


A 


m,m,2 


AfEki 


-5{5 — 27r) 


(33) 


Working similarly with 2 (where it would be necessary to use e < 0 since 2 arises in the evaluation of 

[-^ - ’ ]|5-)^ we would obtain: 


o{d) ^ ^3_ 
AEkd 


5{5 + 27r) 


(p) fp) 

4 m m 1 and 3 ‘^an be evaluated following the same procedure, resulting in: 


aip) ^ oip) 
r^m,m,l r^m,171,3 


J6^ 

87rfc„ 


(34) 


(35) 


Note that in the above mentioned calculations, we have ignored the slope s„ of the surface. Eor a sufficiently fine 
sampling, this approximation is a valid one, since rn- 1/2 ~ i"n+i/ 2 ^ thereby implying « 0. 


3) Eor n = m + 1, the integrand in 


j(p) 


is not singular and can be approximated centroidally: 


p(p) _ ^Iq ikpRm,m+7/4) ^ 
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m+7/4,m 


(36) 
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(v) 

m +1 2 singular and must be evaluated using the procedure described above. It can be shown that in this case, the 
limiting value of the integral is same as the value of the integral at the cylinder surface; 


A 


(p) __ 


Care must be taken in evaluating the limit in /3, 


(p) 


o(p) _ J STTfco 

■^m,m+l,3 I jS 


^(S-2Tr) forp = 0 


+ foi-P = rf 

4) The case of m = n + 1 is identical to the case in which m = n — 1. We give the final results below: 

M ^ forp = 0 

l5&(« + 2>r) for P = d 


o(p) 

1,2 
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(p) 

m,m—1,3 


STrfcd 
Airkp 

Hq ^ {kpRra,m-7/i) 
S k^m,ra—7/4: 


m — 


7/4,r 


m— f/4,m 


Consider the inner product of with Pmi4>), which can be approximated centroidally: 




Pm{4>) d(j) = em^ 5exp^-jk^Tm cos{mS - ^o)) 


(37) 

(38) 

(39a) 

(39b) 

(39c) 

(40) 


Finally consider the inner product of K{(j)) with Pm{4>)- 

271- {m+l/2)S 


^ r 3(5 (5 

KPm{4^) ~ ^ ^ I ~ g(^m—1 “1“ ^m+l) 

n n-0 


(41) 


(m—1/2)5 

Imposing Eq. ([^ by taking the inner product of LHS and RHS with Pm(0) for m G {0,1, 2,... — 1}, we obtain the 

following system of linear equations: 


36 6 

— ~^km H" ^(^m—1 “t“ ^m+l) “t“ ^ ^ ^nmJm 4” ^ ^ 


n—0 

N-1 


n—0 

N-1 


36 6 

0 = + ^(^m—1 H“ ^m+l) + j'Hdkd ^ ^ ^nmjm 4“ ^ ^ Pnm^'n 


(42a) 


(42b) 


n—0 


n—0 


which is a system of 2 A^ linear equations in 2 A^ unknowns kn^ n G { 0 , 1 , 2 , 3... — 1 }) that can be solved to obtain 

jn and kn and thus J((^) and iT(0) through Eq. (j^. Once J(0) and K{(j)) are known, the scattered fields can be evaluated 
through: 

Er = -jmkoAT - ^ ^ ^ 


dx dy J 

wherein and F are defined through Eq. 0- A quantity of interest related to the scatterer is the two-dimensional scattering 
cross-section a{(j)) defined by: 

■ (^.ca )2 


= lim 27rr 


p 2 

^0 


(44) 


Eq being the amplitude of the incident plane wave. The limiting form of as r — > c» can be evaluated using Eqs. 
and Q. Note that as r —>■ oo, 


and thus 


R = i/r^ -I- rs{(t>'Y ~ 2pps(A) cos((/) — (f)') « r — rs{4>') cos{(j) — 4>') 
H^o\koR) « ^J~^^exp{-j{kor - 7 r/ 4 )) exp{jkorsi^') cos{^ - (j)')) 


H^^'^'ikQR) « -j 


nkor 


exp{-j{kor - 7 r/ 4 )) exp(jfcors(())') cos{(j) - (p')) 


(45) 

(46a) 

(46b) 


















Similarly, 


(47) 


F{(j), (j)') « r{rs{(j)') cos{(j) - (j)') + r's{(j)') sin(^' - (f))) 


Using Eqs. 0 and ([43|, we get: 


lim = - 

r—¥oo 


2 fcp 
Trfcgr 4 


exp(-jfco(r - 7r/4)) 




- K{4)'){rs{(t)') cos{(j) - (j)') + r's{(t)') sin(0' - 4))) exp(jfco?'s(<ii’') cos((/> - cj)') d(j)' 


and thus 




+ - A'(^')(rs(0')cos(0 - 4 >') + rs(i^') sin(0' - ^)) 


expO’/co^'s(<?!>') cos((/i - (/>') d(i>' 


(48) 


(49) 




Fig. 3: Plot of a{4>) computed using MoM and compared with the analytical solution for different discretizations of the cylinder 
surface. In all the calculations, a = 2Ao, permittivity = 2eo, and permeability fici = /ig- 


Using Eq. 0, a{cj)) can be expressed in terms of j„ and kn- 
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where 


an (</>)=??0 / rl{(j)’) + r|'(</)') exp(jA:o?’s(<?!>') cos((/) - 0')) 


(51a) 


&„((?!)) = / A„(( 5 ()')(»'S(</>') cos(0 - (/>') + r's{(j)') sin((;i' - (j))) exp(j/co^s(</>') cos((/i - 0')) 


(51b) 


Since the integrals in Eq. are not singular, they can be approximated centroidally; 

5 n ,,, ,, „,. ,, 3^ 


ani4>) ~ Vo 


+ 4-1 exp(jfcor„- 3/4 cos((n - 3/4^) - 4>) + -^y'rl +six 


exp(jfcor„ cos{n8 - 0) + ^^^’^+ 3/4 + s^+i exp(jfcor„+ 3/4 cos((n + 3/45) - 0) 


(52a) 


52 

~ ^(’’"- 3/4 cos((n - 3/4)5 - </) + s„_i sin((n - 3/4)5 - 0))x 

O 

3^2 

exp(jfcor’„_ 3/4 cos((n — 3/4)5 — (/)) H—cos(n5 — (/) + s„ sin(n5 — </)) exp(jfcorn cos(n5 — (/)) 

52 

+ -^(r„+ 3 / 4 C 0 s((n + 3/4)5 - </) + s„+i sin((n + 3/4)5 - ^)) exp(jA:o?'n+ 3/4 cos((n + 3/4)5 - </)) 


(52b) 


Fig.|3] shows a comparison between the MoM results and analytical calculation based on the Mie theory Ch. 11.5] for 
a smooth cylinder rs{(j)) = a = const. It can be seen that for 5 ~ the MoM described above can be used to accurately 

compute the scattered fields. 


III. Scattering from a perfectly conducting cylinder with axial roughness 

Consider a perfectly conducting cylinder with it’s axis oriented along the 2 : axis and the surface being described by r{(j), z) = 
rs(z) = a + h(z), where h{z) is the roughness function and a is the mean radius of the cylinder. We consider an azimuthally 
symmetric incident field given by: 

E“^(r,(/,z) =£;™^(r,z)(/ (53) 

The above field corresponds to an incoming cylindrical wave incident on the cylinder. To compute the scattered fields, we use 
the EFIE at the surface of the PEC Q Ch. 1.7]: 


nxE* 


= X (VV • A + k^A) 
ko 


where 


A(r) = J(r') 


,^exp(-jA:o|r-r'|)^3^, 


47 r|r — r'l 


(54) 


(55) 


Symmetry dictates J(r) = J{z)(f>. Let r = (r, (j), z) and r' = (r', (/', z'). Noting that J is actually a surface current, Eq. 
be written as (R = [r^ + rs(z')^ — 2rrs{z’) cos{(j) — (j)') + {z — z')^]^/^): 


55 


00 27 r 


A(r) = J J rs{z')^Jl + r'^{z')d(j)'dz' 


(56) 


2'=—00 d )'—0 


Observing that </' = — sin (j)' x + cos (/' y and performing the substitution 0 = (/' — </, we obtain (since the integrand is a 
periodic function of 9 or </', the exact limits are irrelevant and integration over any period will give the same result, thus the 
limits have been chosen from — tt to tt for convenience.): 


00 TT 

A(r) = J dz'J{z')rs{z')^Jl + r'g{z') J {—sm{(j3 + 9)x + cos{(j) + 9)y) 


exp(-jfcoi 7 ) 

47ri? 


dO 


(57) 


Also note that R is an even function of 9. Consider the integral over 9: 


(- sin(0 + 9)x + cos((^ + 9)y)d9 = ^ [ cos9 ^''^[ d9 (58) 


47ri? 


47ri? 


9 — — 7\ 
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with (j) = — sin^ x + cos </> y. To arrive at this result, we expanded sin((/) + 6) 
cos (j) cos 6 — sin (p sin 9 and used the fact that; 


Sin o'- dO = 0 

AttR 


9 — — 7T 


sin (p cos 9 + cos p sin 9 and cos{p + 9) = 


(59) 


since the integrand is an odd function of 
the integral, we obtain: 


Plugging Eq. (581 into Eq. (571, and noting that is a constant with respect to 


A{r) = A{r, z)(p = (p / dz'J{z')rs{z')\/1 + rg{z') / cos9 


^ex p{-jkoR) 
AttR 


d9 


z'——oo 


9 = — 7\ 


Note that Eq. (601 immediately implies that V • A = 0. Therefore, the EEIE (Eq. ([54|) simplifies to; 


z) 


= 3 mkoA{r, z) 


(60) 


( 61 ) 


To numerically solve Eq. (611, we consider only a hnite length of the cylinder, say from z = —L/2 to z = L/2. As in Section 
[H] a piecewise linear approximation to the cylinder surface between N + 1 discrete points is assumed. Let the points be labeled 
as {zn, r„), where n G {0,1,2,... TV} and Zn = nS, r„ = rs{zn), with 6 being the discretisation length S = L/N. The radial 
coordinate of a point on the cylinder surface at z, where Zn < z < Zn+i is approximated by: 


/ \ -^n+l Zn Zn^l Z'^ 

rs[z) « --+ 


■ = Zn+i/2 + s„ + 0 6^ (62) 


where we have dehned z 


n+l/2 


2 ' S 

z„+i+z„ _ z^^^-zr, ^ expand the unknown J(z) on the pulse basis functions: 


N-l 


J{z)^ jnPn+l/2{z) 


(63) 


n=0 


and test Eq. (611 through the delta testing functions T^+i/^iz) = S(z — Zn+i/ 2 )- Consider, now the inner product of each term 
in Eq. (611 with T„+i/2(z). Eor instance: 


E;“=(r,z) 


Pm+l/ 2 (z) dz — ^m+1/2) — Gn 


(64) 


where rm+ 1/2 = rs{Zm+i/ 2 ) « Similarly: 


n 00 TT 

J A{r,z) dz = J dz' J{z')rs{z')^ I + r'^{z') J cos9 

0 2'= —00 9— — 7T 


exp(-}fcoi? m ) 


AirRr, 


N-l 

dO — ^ ^ C^m,njn 
n—0 


where 


^m.n — 


00 TT 

J dz'Pn+i/ 2 {z')rs{z')^l + r'^{z') J cos9 


exp(-jfcoi? m) 


47ri?„ 


d9 


Z— — 00 

(n+l)S 


J dz'rs{z')\Jl + r'g{z') J cos 9 

'—n6 9— — '7 t 


9 = — '7T 

exp(-jfcoJ7 

m) 

AirRjn 


d9 


(65) 


( 66 ) 


with Rjn = R = [rs(z')^ + ?’m+i /2 “ 2rs(z')r„+i/2 cos(^ - </>') + (z' - Zm+ 1 / 2 )^]^/^. We treat the evaluation of am,n 
seperately for two cases; 

1) Eor m ^ n, the integral in z' can be approximated centroidally resulting in; 

TT 

X fT ,—T f aexp(-jfco-Rm.n) 

am.n = d?'n+l/2\/l + Sn / COS ^-——- d9 (67) 

9=—tt 


where = [rYi /2 + ’'n+1/2 + - 2r„_,_i/2fm+i/2 008 0]^/^. The integral over 9 in a.m,n would have to 

be evaluated numerically using quadrature methods. 
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2) For m = n, the integral is singular and the singularity would have to be integrated out. Following the method described 
in ij^ Ch. 8.7], we perform a change of variables z' ^ z' — {m 1/2)6 and write am, 7 n 

6/2 -TT 

exp{-jkoR n 


j dz'rs{z' + (m + l/2)5)^Jl + rg[z' + (m + l/2)<5) J cos9 

9 — — 7T 

TT 

= 4 f dz'rs{z' + {m + 1/2)8)^ 1 + r'g{z' + (m + 1/2)5) f 


z' — — 6/2 


47ri?„ 


AirKn 


dO 


de 


z'=0 


e=o 


(5/2 TT 

4r„+i/2\/l + J dz' J cos 9 


= 4r„+i/2\/l + 4 


exp(-jfcoi? 

rj 

471-7?^ 

<5/2 TT 

f f exp{-jkoRjn) - 1 , , 


d6» 


s /2 


R-n 


dz d9 


R„ 


-d9 


z'=o e=o 


z'=o e=o 


1 


In 


^+l/2 


8r„ 


(+1/2 


S/2 


dz' 


1 


( 1 + 1/2 


In 


8r„ 


(+1/2 


dz' 


( 68 ) 


z'=0 

The first integral is not singular and can be evaluated numerically. The second integral is also not singular can be rewritten 
as: 

s /2 

( r ^ 1 / z' N t 

dz' 


z'=o e=o 
s /2 


—— d9 -In 

Ttm '^m+1/2 


8 r„ 


1 + 1/2 


/ 


2 dz' 


4r£,, 1 /T + z'^ 


E 


\ 


4j.2 


+ 1/2 


4^.2 


1 + 1/2 


+ z ■ 


0 V "1+1/2 

where E{x) is the complete elliptic integral of the first kind, defined by: 

ir/2 


( 1 + 1/2 


In 


8 r„ 


1 + 1/2 


(69) 




d/3 


— X sin /3 


(70) 



Fig. 4: Plot of |i3^“(r, z)| as a function of z computed using MoM for an incident field of the form 73™'^ = H/\kQr). In 
this simulation a = 2Ao, r = 2.5Ao, S = O.IAq and L = ISAq. The analytically computed fields for an infinitely long cylinder 
for the same dimensions is also shown to illustrate the diffraction effects due to the finite length of the cylinder. 
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The integral in Eq. 


can be evaluated numerically. The last integral can be analytically evaluated to yield; 

(5/2 


1 


z'=0 


^7n+l/2 


■In 


8 r„ 


(+ 1/2 


dz' = 


2r„ 


1+1/2 


In 


16r„ 


1+1/2 


- 1 


Finally imposing Eq. ([M), we obtain; 


N-l 

— jVO^O ^ ' Hm.rijn 
n=0 


(71) 


(72) 


which can be solved to obtain jn- Once jn are known, the scattered fields at any point (r, z) can be evaluated. For instance, 
E‘^^°‘{r,z) is given by; 

z) = -jko7joA{r, z) 


where 


OO TT 

A{r,z)= j dz'J{z')rs{z')^ 1 + r'g{z') J cos9 


exp{-jkoR) 

inR 


dO 


z' = —OO 




Substituting for J(z') from Eq. (63 i, we obtain; 


N-\ 


= -jkoTjo 


n—0 


where 


(n+l)(5 


J dz'rs{z')\Jl + rg{z') J cos 9 

'—n5 6— — ‘7 t 


Sr, 


n+l/2 


x/T" 


COS0 


exp{-jkoR) 

AttR 

exp{-jkoR„ 

ATvRn 


d9 


d9 


(73) 


(74) 


(75) 


(76) 


with Rn = Ayi’^+i /2 + - 2rr„+i/2 cos0 + (z - (n + 1/2)5^. 

It is not possible to simulate the response of an infinitely long cylinder with axial roughness, since the MoM would require 
us to truncate the cylinder in the z direction. Truncation of the cylinder results in diffraction effects from the cylinder edges 
(as is shown in Fig. |^, which is undesirable if the objective of this exercise is to simulate the response of an infinitely long 
cylinder. However, a standard approach to this problem is to use a tapered incident field Q. For instance, consider an incident 
field of the form; 


jpinc _ 

Es - Wo 


exp - 


k'^Wn 


H^\'r^k‘l — fc2)exp(jfcz) dk 


(77) 


Note that H\^\r^kQ — exp(jfcz) is a cylindrical wave propagating from infinity to the axis of the cylinder and with 
a component kz of the wave-vector along the z direction. By adding such waves with a gaussian envelope exp(—A:^Wg/4), 
we have constructed an incident field which is approximately gaussian in the z coordinate. This can be explicitly seen if 
koWo » 1, in which case only very small k contribute significantly to the integral in Eq. (77 1 . This allows us to make the 
following Taylor series expansion to the second order in k: 

H\^\r\lkl-k^) - H\^\kor) - (kor) 


(78) 


The integral in Eq. (|77|) can then be evaluated analytically to yield; 


= 2v^exp - 


Wr 






(1)' 


(kor) 


(79) 


Thus, if L ^ wq, then the cylinder behave almost as if it is an infinitely long cylinder as far as the gaussian incident field 
is concerned. Fig. shows the comparison between analytical results and MoM results for a smooth cylinder for iV = 50 and 
N = 150. Clearly, <5 ^ ^ is sufficient to accurately simulate the response of the cylindrical scatterer. 
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Fig. 5; Plot of z)\ as a function of 2 computed using MoM and compared with the analytical solution for different 

discretizations of the cylinder surface. In all the calculations, a = 2Ao,r = 2.5Ao, wq = 5Ao and L = ISAq. 
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